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Abstract. Recent VLT SINFONI observations of the close environments (~ 30 pc) of nearby AGNs have shown that thick 
gas tori and starbursts with ages between 10 and 150 Myr are frequently found. By applying these observations to a previously 
established analytical model of clumpy accretion disks, we suggest an evolutionary sequence for starburst and AGN phases. 
Whereas the observed properties of the gas tell us about the current state of the torus, the starburst characteristics provide 
information on the history of the torus. In the suggested evolution, a torus passes through 3 different phases predetermined by 
an external mass accretion rate. Started by an initial, short, and massive gas infall, a turbulent and stellar wind-driven Q ~ 1 
disk is formed in which the starburst proceeds. Once the supernovae explode the intercloud medium is removed, leaving a 
massive, geometrically thick, collisional disk with a decreasing, but still high-mass accretion rate. When the mass accretion rate 
has significantly decreased, the collisional torus becomes thin and transparent as the circumnuclear disk in the Galactic center 
of the Milky Way. Variations on this scenario are possible either when there is a second short and massive gas infall, in which 
case the torus may switch back into the starburst mode, or when there is no initial short massive gas infall. All observed tori up 
to now have been collisional and thick. The observations show that this phase can last more than 100 Myr. During this phase 
the decrease in the mass accretion rate within the torus is slow (a factor of 4 within 150 Myr). The collisional tori also form 
stars, but with an efficiency of about 10 % when compared to a turbulent disk. 

Ksy words. Galaxies: active - Galaxies: nuclei - ISM: clouds - ISM: structure - ISM: kinematics and dynamics 



1. Introduction 

In the unification scheme for active galactic nuclei (AGN) the 
central massive black hole is surrounded by a geometrical thick 
gas and dust torus (see, e.g., Antonucci 1993). If the observer's 
line-of-sight crosses the torus material, the AGN is entirely ob- 
scured from near-IR to soft X-rays and only visible at X-ray en- 
ergies if the gas column density is not too high (Seyfert 2 galax- 
ies). On the other hand, if the torus is oriented face-on with 
respect to the observer, the central engine is visible (Seyfert 
1 galaxies). The spectral energy distributions (SEDs) of most 
quasars and AGN in Seyfert galaxies have a pronounced sec- 
ondary peak in the mid-infrared (mid-IR) (e.g. Sanders et al. 
1989; Elvis et al. 1994), which is interpreted as thermal emis- 
sion by hot dust in the torus. The dust is heated by the primary 
optical/ultraviolet (UV) continuum radiation, and the torus ex- 
tends from the dust sublimation radius outwards (Barvainis 
1987). 

The geometrical thickness of the torus in the gravitational 
potential of the galactic nucleus implies a vertical velocity dis- 
persion of about 50-100 kms -1 . If one assumes that the disk 
is continuous, i.e. thermally supported, this corresponds to a 
temperature of ~ 10 5 K. Since this is beyond the dust subli- 
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mation temperature (~ 10 3 K), thick tori have to be clumpy 
or must be supported by additional forces other than thermal 
pressure. Krolik & Begelman (1988) proposed a clumpy torus 
model where the clumps have supersonic velocities. Vollmer et 
al. (2004) and Beckert & Duschl (2004) elaborated this model 
in which orbital motion can be randomized if magnetic fields 
permit the cloud collisions to be sufficiently elastic. Vollmer 
et al. (2004) found that the circumnuclear disk (CND) in the 
Galactic center (Giisten et al. 1987) and obscuring tori share 
the same gas physics, where the mass of clouds is in the range 
20 - 50 Mq and their density close to the limit of disruption 
by tidal shear. A change in matter supply and the dissipation of 
kinetic energy can turn a torus into a CND-like structure and 
vice versa. Any massive torus will naturally lead to sufficiently 
high mass accretion rates to feed a luminous AGN. 

If and how efficient these clumpy tori form stars is an open 
question. The large majority of observational studies probed 
the nuclear star formation on scales of a few hundred parsecs 
(see, e.g., Sarzi et al. 2007, Asari et al. 2007, Gonzales Delgado 
& Cid Fernandes 2005, Cid Fernandes et al. 2004). These stud- 
ies resulted in a general view that about 30 % - 50 % of the 
sample AGNs are associated with recent (ages less than a few 
100 Myr) star formation on these scales. Thanks to the high 
spatial resolution of the near-infrared adaptive optics integral 
field spectrograph SINFONI, it has only recently become pos- 
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sible to study the environments of AGN on the 10 pc scale. 
Davies et al. (2007) analyzed star formation in the nuclei of 
nine Seyfert galaxies at spatial resolutions down to 0.085". 
They found recent, but no longer active, starbursts in the central 
regions which occurred 10 - 300 Myr ago. Moreover, Hicks et 
al. (2008) were able to measure the rotation and dispersion ve- 
locity of the molecular gas in these galaxies using the 2.12 
H2 (1-0)S(1) line. Surprisingly, all molecular gas tori have high 
velocity dispersions and are therefore geometrically thick. 

In this article we compare the observations of Davies et al. 
(2007) and Hicks et al. (2008) with the expectations of ana- 
lytical models of clumpy accretion disks developed in Vollmer 
& Beckert (2002, 2003) and Vollmer et al. (2004). In a first 
step, we test whether these models are able to describe obser- 
vations. In a second step, these models allow us to investigate 
the scenario hypothesized by Davies et al. (2007) where spo- 
radic, short-lived starbursts are due to short massive accretion 
events in the central region, followed by more quiescent phases 
until there is another episode of accretion. This scenario is cor- 
roborated by a closer look at the nucleus of NGC 3227 where 
Davies et al. (2006) found signs of a past starburst (~ 40 Myr 
ago) and a presently quiescent gas torus with a Toomre param- 
eter Q > 1. 

2. The theory of clumpy gas disks 

Within the framework of Vollmer & Beckert (2002, 2003) and 
Vollmer et al. (2004) clumpy accretion disks are divided into 
two categories: (i) turbulent and (ii) collisional disks. In case 
(i) the ISM is regarded as a single entity which changes phase 
(molecular, atomic, ionized) according to internal (gas density, 
pressure, magnetic field) and external (gravitation, radiation 
field, winds) conditions. Energy is injected into a turbulent cas- 
cade at the driving length scale (large scale) and dissipated at 
the dissipation length scale (small scale). We identify the dis- 
sipation length scale with the characteristic size of selfgravitat- 
ing clouds. These clouds decouple from the the turbulent cas- 
cade and constitute the first energy sink. The source of energy 
which is injected at the driving scale to maintain turbulence can 
be either (i) mass accretion in the gravitational potential of the 
galactic center (fully gravitational FG model) or (ii) supernova 
explosions (SN model). In the collisional case energy is also 
supplied in the process of mass accretion in the gravitational 
potential of the galactic center and dissipated via partially in- 
elastic cloud-cloud collisions. The actual dissipation rate in in- 
dividual collisions is largely unknown. The disk evolution is 
mainly driven by the external mass accretion rate. Since these 
models are equilibrium models, we assume that the mass ac- 
cretion rate is constant throughout the region of interest when 
averaged for a sufficiently long time (~ here 51 is the 

angular velocity of circular orbits in the gravitational potential 
of the galatic nucleus) and maintained for at least the turnover 
timescale i?/t>t U rb °f S as m me disk, where R is the distance 
from the center of the galaxy and uturb the characteristic speed 
of turbulent eddies. All models give access to the global param- 
eters of the disk and the local parameters of the most massive 
clouds (see TableHJ. The free parameters of the models are the 
Toomre parameter Q, the disk transparency $ (Eq. |7), and the 



Table 1. Model parameters and their meaning 



large scale 


disk 


R 


galactic radius 




rotation velocity 




angular velocity 


M dyn 


total enclosed (dynamical) mass 


M gas 


total gas mass 


^turb 

Q 


gas turbulent velocity dispersion 
Toomre parameter 


p 


midplane density 


E 


surface density 


H 


disk height 


M 


disk mass accretion rate 


V 


viscosity 


'driv 


turbulent driving length scale 


s 


scaling parameter between driving 
scale length and cloud size 




cloud volume filling factor 


c 


viscosity scaling parameter 


,H 

% 


disk vertical free fall time 


p* 


star formation rate per unit volume 


V 


hLdl 1UI IIldjLKJIl Idle JJcI UI11L bLlildjUJ 




star formation rate 


£ 


conversion factor for SN energy flux 


T] 


star formation efficiency in the 
collisional model 


i 


linking factor between Mgas and M 
for torus evolution 


■3 


disk transparency 




mass accretion rate onto the black hole 


Jvl wina 


wind outflow rate 




AGN UV radiation field 




ClOLlLlb 


icon 


timescale for cloud — cloud collisions 


/coll 


cloud mean free path 


M cl 


cloud mass 




cloud radius 


N d 


cloud surface density 


C s 


local sound speed within the clouds 


ts 


sound crossing time of clouds 


,cl 
Iff 


cloud free fall time 


T 


temperature 


Ci 


sound crossing time of the ionized gas 



mass accretion rate M . Other parameters are fixed using the 
Galactic values (Vollmer & Beckert 2002, 2003). Each model 
has an associated star formation rate. In the following we de- 
scribe these models in more detail. 

2.1. Turbulent disks 

In Vollmer & Beckert 2002 (Paper I) we developed an analyti- 
cal model for clumpy accretion disks and included a simplified 
description of turbulence in the disk. In contrast to classical ac- 
cretion disk theory (see, e.g., Pringle 1981), we eliminated the 
"thermostat" mechanism, which implies a direct coupling be- 
tween the heat produced by viscous friction and the viscosity 
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itself. The viscosity is usually assumed to be proportional to 
the thermal sound speed. Thus, the (gas) heating rate depends 
itself on the gas temperature. This leads to an equilibrium cor- 
responding to a thermostat mechanism. Instead, we use energy 
flux conservation, where the potential energy that is gained 
through mass accretion and differential rotation is cascaded by 
turbulence from large to small scales and dissipated there. 

One fundamental approximation is that the kinetic energy 
is dissipated (removed from the turbulent cascade) when the 
gas clouds become self-gravitating. Turbulence transfers the 
energy from the driving wavelength i<jriv to the dissipation 
wavelength Z^iss, which corresponds to the size of the largest 
selfgravitating clouds. The two length scales are linked by the 
scaling parameter f . For a Kolmogorov-like turbulent energy 
spectrum £ = (/driv/'diss) 3 - In addition, the modeled disks 
have a constant Toomre-Q parameter: 



Wturb ^ 

ttGS 



> 1 



(1) 



where i>t U rb is the turbulent velocity, ft the angular velocity, G 
the gravitational constant, and £ the gas surface density of the 
disk. If we can approximate the total gas mass within a radius 
R by M gas = 7ri? 2 S the Toomre parameter can be rewritten 



Wturb Afdyn 
^rot.K M„ as 



(2) 



where Afdyn is the total enclosed mass and w rot the Keplerian 
rotation velocity. 

Furthermore, we use the following prescription for the vis- 
cosity: 



(3) 



We obtained a set of equations with 3 free parameters: the 
Toomre parameter Q, the scaling parameter £ > 1, and the 
mass accretion rate within the disk M. The mass accretion rate 
is the gas mass transported per time through the gas disk from 
large to small radii. It must be supplied from the galaxy at the 
outer radius of the turbulent gas disk and is constant at all radii 
in the disk. This set of equations can be solved analytically and 
the results were already used in Paper I to describe properties of 
our Galaxy. The solutions depend on the parameters Q, £, M, 
v Iot , and R. It turned out that the driving wavelength equals the 
disk height Zdriv = H. 

In a second step we included the energy input due to su- 
pernova (SN) explosions (Vollmer & Beckert 2003, Paper II). 
The energy flux provided by SNe is transfered by turbulence 
to smaller scales where it is again dissipated. The SN energy 
flux is assumed to be proportional to the local star formation 
rate. The local star formation rate is taken to be proportional 
to the mean density and inversely proportional to the local free 
fall time of the clouds. These clouds have sizes that are a factor 
8 smaller than the driving length scale. The factor of propor- 
tionality is the probability to find a self-gravitating cloud, i.e. 
the volume filling factor. The integration length in the vertical 
z direction is assumed to be the turbulent driving scale length, 
i.e. the vertical height in the disk where self-gravitating clouds 
are found: £* = Zdriv The SN energy per unit time _Esn per 



area AA is therefore proportional to the local star formation 
rate £*: 



Esn 
~~AA 



= ££, 



(4) 



where the factor of proportionality £ is independent of the ra- 
dius in the disk. Its normalization with Galactic values yields 
f = 4.6 10" 8 (pc/yr) 2 . 

In the FG model the energy transported through the turbu- 
lent cascade is supplied by mass accretion, which leads to an 
energy flux equation of the form 



pv 



7 turb 



Idxi 



i • on 

MV Tn t — — 

2vr dR 



(5) 



In the case of SN driven turbulence the energy flux is deter- 
mined by the star formation rate 



pv 



y turb 



(6) 



Furthermore, we take into account that the clouds are formed 
during the interaction between SN remnants at the compressed 
edges. The size of the clouds l c \ is consequently smaller than 
the turbulent driving wavelength and we use 6 = idriv/^cl with 
6>l. 



2.2. Collisional disks and the stability of the clouds 

If the selfgravitating clouds are stable, their collisions will give 
rise to angular momentum redistribution again described by 
an effective viscosity. An equilibrium disk can be formed if 
there are fragmenting collisions or partially elastic collisions 
(the clouds are supposed to be magnetized). If the collisional 
timescale t co \y is longer or equal to the dynamical timescale, 
the resulting viscosity can be written as 



1? Wt ur b-ff 



(7) 



where the disk transparency d = t co i\Cl > 1 and H is the disk 
height. It follows that the collisional energy dissipation rate is 



AE 
AAAt 



f 



^turb 



^turb 



coll 



/ 



coll 



■OH 



(8) 



Here the factor /, which has been omitted in Vollmer et al. 
(2004), accounts for the mean fraction of cloud mass partici- 
pating in the highly supersonic cloud collisions. For constant 
density clouds Krolik & Begelman (1988) argue that / = 0.2. 
For more centrally condensed, self-gravitating clouds we will 
use a factor / = 0.1 in this paper. This geometric factor / 
not only reduces the average energy dissipation rate in cloud 
collisions but also the angular momentum redistribution in the 
collisions. The reduced energy dissipation is accompanied by 
a correspondingly reduced mass accretion rate M for the same 
transparency in the model disks v — > fv in Eq.((7]l. 

Since in the FG model i<jriv = H, the FG model and the 
collisional model are formally equivalent. We can thus use 
Eq. ([TJ-O replacing £ by ■d to describe the collisional model. 
Nonetheless the interpretation of £ and d is completely differ- 
ent. 
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The cloud size r c \ and the volume filling factor $y of 
clouds can be derived using their mean free path (see Vollmer 
et al. 2004) 



the characteristic timescale for the cloud collapse, i.e. the non- 
averaged local free fall time tg m . 



Zcon =#H = 



4r, 



cl 



3$x 



and the fact that the clouds are selfgravitating 



.cl 



/37T$V 

32G> 



= t„ = 



(9) 



(10) 



P* « Id 



(15) 



Since tg oc p » this corresponds to a Schmidt law of the form 
p* oc p3 . The factor of proportionality is given by the probabil- 
ity to find a self-gravitating cloud, i.e. the volume filling factor 
4>y. Thus, the star formation rate is given by 



where t g l is the free fall timescale within clouds, p the disk 
overall gas density, t s the sound crossing timescale, and c s the 
sound speed. This leads to 



(16) 



r c i = 
and 



Qcl 



8 Slvturb^ 



Qcl 



6 <4 rb ^ 2 



(11) 



(12) 



The cloud mass is then 



M c i = 



v P r d = 



64 



sG^Qs/^-sr 1 .(13) 



These clouds move supersonically within the intercloud 
gas. If the intercloud gas is ionized, typical Mach numbers are 
about 10. Therefore, we expect that the clouds might be de- 
stroyed by Rayleigh-Taylor instabilities. However, Rodiger & 
Hensler (2008) showed that these instabilities are suppressed in 
the presence of (i) a sufficiently strong gravitational field of the 
clouds or (ii) a strong magnetic field. The ram pressure exerted M* = r]M sas 'd~ 1 Q, . 
on the clouds is p ram = pintv^, wnere Pint is the intercloud 
gas density and v c \ is the cloud velocity. Following Rodiger & 
Hensler (2008) Rayleigh-Taylor instabilities are suppressed if 
the gravitational acceleration g is higher than the acceleration 
due to the drag by ram pressure an- This yields approximately: 



Furthermore, we assume that stars are only born in the mid- 
plane of the disk in regions that have the size of the turbulent 
driving length scale Zdriv, because the clouds can collapse only 
within the turbulent timescale i turb = ldmv/v tm b- We thus ob- 
tain 

S„ = p» ^driv (17) 

for the mass surface density turned into stars. 

3.2. Collisional disks 

For the collisional disk we assume that the star formation rate 
is proportional to the overall density p and the cloud collision 
frequency t~^ n = i9 _1 fh 



(18) 



where i] is an a priori unknown efficiency factor. In terms of 
stellar mass per time this gives 



(19) 



4. Thick disks in a generic galactic center 



M C \G fern 
r r\ ^Cl 



(14) 



For the intercloud gas density we use the value of the Galactic 
Center p int = 10 3 cm~ 3 (Erickson et al. 1994). Using M c \ = 
10 M Q , r c i = 0.02 pc, v c \ = 100 kms -1 , the gravitational ac- 
celeration is about 4 times higher than the ram pressure drag. 
If the intercloud gas has a 10 times higher density, only a mag- 
netic field with a field strength of B ~ ylO Pint^i ~ 1 m ^ 
can stabilize the clouds. This kind of field strength is observed 
in the Circumnuclear Disk in the Galactic Center (Plante et al. 
1995). It is thus plausible that the torus clouds are stable against 
Rayleigh-Taylor instabilities. 

3. Star formation in clumpy gas disks 

3.1. Turbulent disks 

Following Paper II we assume that the star formation rate is 
proportional to the mean density of the disk and the inverse of 



The theoretical model described in the previous sections 
can now be applied to recent near-IR high-spatial resolution 
SINFONI observations of nearby AGN. Davies et al. (2007) 
showed that there had been recent star formation in the cen- 
tral few tens of parsecs in a sample of nearby AGN. Following 
on from this, Hicks et al. (2008) showed that distribution and 
kinematics of the central concentrations of gas were simi- 
lar to those of the stars, and was geometrically thick with 
Wturb/^rot > 1/2. They argued that this gas comprised the 
large scale structure of the tori, implying that tori can form 
stars. And they gave an estimate for the gas mass in this region 
as 10% of the dynamical mass. Before we derive the physical 
disk parameters for the individual AGNs observed by Davies 
et al. (2007), we give an overview over the different types of 
thick tori. Note that z9 is the disk transparency, i.e. $ < 1 im- 
plies an opaque disk whereas a large i? results in a transpar- 
ent disk. For a generic galactic center we assume a dynamical 
mass (which for radii R larger than a few pc is dominated by 
the stellar content), of Md yn = 10 8 M Q and a gas mass of 
M gas = 10 7 M©, both within a radius of R = 20 pc. The en- 
closed mass leads to a rotation velocity of 150 kms -1 at this 
radius. In addition, we adopt a cloud internal sound speed of 
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c s = 1.5 km s _1 . The value for the sound speed corresponds to 
a cloud temperature of ~ 500 K when only thermal gas pres- 
sure is considered. The sound speed is a measure of the pres- 
sure support against self-gravity and additional contributions 
to the pressure gradient inside clouds like magnetic fields may 
contribute. The adopted sound speed leads to cloud masses of 
M c \ oc c 4 ~ 10 Mq. The disk is stable with a Toomre param- 
eter of Q = 4.7 for a turbulent velocity supporting the vertical 
thickness of Ut ur b = 70km/s. Moreover, we assume a star for- 
mation rate of M = 0.1 M yr _1 . This can be considered an 
upper limit to the current star formation rate in the central few 
tens of parsecs. Davies et al. (2007) showed that while star for- 
mation had occured there recently, it has now ceased. Based on 
the exponentially decaying starburst model they used and the 
time averaged star formation rates they estimated, the current 
rates are expected to be below this limit. 

In the following we generate a generic set of models which 
(i) can describe different evolutionary phases of a ~ 20 pc 
scale gas disk in terms of mass accretion rate, disk thick- 
ness/transparency, and star formation rate and (ii) can repro- 
duce the observations at a generic time in the evolutionary se- 
quence. We first compare different disk models, representing 
evolutionary stages, to the observations of the current state of 
the disk (Sec. 14. 1144. 5b . It is shown that a massive accretion 
event leads to a relatively thin turbulent disk which forms stars 
at a high rate (SN model). Once the SN explode the intercloud 
medium is blown out, leaving only dense, compact clouds. This 
results in a collisional disk that may assume one of several 
states (Sec. 14.2144.51 ) depending on the external mass accre- 
tion rate M, the Toomre Q, and transparency -d parameters. In 
Sec. [5] we then investigate the evolution of the collisional disk. 

4.1. A massive turbulent disk creates a starburst 

We first try applying the turbulent SN disk model. With the 
parameters described above it would have a mass accretion rate 

M = 2.18 A/-X 2 as«t 4 urb^ 2 r 1 (20) 

in excess of 10 3 M Q yr~ 1 and a ratio between the turbulent 
driving and dissipation length scales of 

(5=1.17 AC 3 G 3 Af g 6 as x; turb ir 6 r 2 > 10 4 . (21) 

Such a model can be discarded, because its mass accretion rate 
is far too high with respect to any supply from the outer galaxy 
and the lifetime M gas /M — 10 4 yr is shorter than the dynam- 
ical time f2 — 1 ~ 10 5 yr. 

On the other hand, if we assume that a turbulent SN disk is 
responsible for the starburst with a star formation rate of = 
1 MQyr -1 (comparable to the initial rate inferred by Davies 
et al. 2007) and that the turbulent velocity was lower («turb = 
20 kms^ 1 ) than the present collisional disk, we find a mass 
accretion rate of M = 2 Moyr -1 and S = 23. These values are 
close to the values for the large-scale Galactic disk (see Paper 
II). A typical cloud has a radius of r c \ — 0.1 pc and a mass of 
M c \ = 10 4 Mq. This kind of disk contains about a thousand 
clouds within R = 20 pc. 

We conclude that a viable turbulent massive gas disk has a 
rather low turbulent velocity and is therefore moderately thin 



with H/R = fturb/^rot ~ 0.13. It then yields large star for- 
mation rates of the order of one solar mass per year. This repre- 
sents a starburst which subsequently will destroy the disk once 
the supernovae explode after about 10 Myr. These explosions 
do not cause any harm to the densest and most massive clouds, 
but they clear the space between the clouds, i.e. they remove the 
initial intercloud medium. We are then left with a collisional 
disk. 

4.2. A massive, opaque, collisional disk 

(Q = 5, = l) 

The collisional disks can be distinguish by their Q and d pa- 
rameters. We start with collisional and opaque disks, i.e. the 
mean free path of the clouds is about the height of the disk H. 
This implies that along a vertical path through the disk there 
is on average one intervening cloud. Along a path in the disk 
midplane towards the center there are typically N ~ 10 clouds 
blocking the direct view. Collisions are frequent in such a torus 
or disk. This yields a mass accretion rate of 

Af = 2^| = 3.4Af Q yr- 1 . (22) 

The volume filling factor of the clouds is $v = 0.004, the 
clouds have typical radii of r c \ = 0.02 pc, and masses of M c \ = 
15 M . 

4.3. A massive, transparent, collisional disk 

(Q = 5,ti = 10) 

The large mass accretion rate in the above model (Sec. 14.21 ) de- 
pends linearly on the collision rate and is a consequence of the 
low transparency. If the disk is more transparent, d — 10, the 
mass accretion rate is therefore M = 0.3 M Q yr~ 4 , the volume 
filling factor of the clouds is $v = 4 10~ 5 , and the clouds have 
typical radii of r c i = 0.002 pc and masses of M c \ = 1.5 M Q . 
We see that for the same Toomre-Q a larger transparency im- 
plies smaller and less massive clouds with smaller volume fill- 
ing factors. 

4.4. A light, opaque, collisional disk 

(Q = 50,$= 1) 

If clouds are large and less dense we can have a disk which is 
light but still optically opaque due to dust in the clouds. For this 
disk class we assume a gas mass of only 1% of the dynamical 
mass, i.e. M gas = 10 6 Mq. This yields a mass accretion rate 
of M = 0.3 Moyr" 1 . The volume filling factor is $ v = 0.04, 
and the typical cloud radii and masses are r c \ = 0.2 pc and 
Mci = 150 M Q . 

4.5. A light, transparent, collisional disk 

(Q = 50, ■& = 10) 

The last disk class is transparent, $ = 10, and has a small gas 
mass, M gas = 10 6 M Q . It has the lowest mass accretion rate, 
M = 0.03 M Q yr _1 , and a low volume filling factor, $v = 
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4 10~ 4 . The cloud radii and masses are the same as those of the 
massive, opaque disk (r c \ = 0.02 pc and M c y = 15 M Q ). 

We conclude that (i) a high mass accretion leads to a mas- 
sive, opaque disk and (ii) a (Q = 5, $ — l)-disk and a (Q = 50, 
{) = 10)-disk share clouds of the same mass and size. These 
clouds are very similar to those found in the Galactic Center 
(see Vollmer et al. 2004). We thus can draw up a disk evolution 
in which a massive, opaque disk evolves with time into a light, 
transparent disk. 

For the rest of the article we assume that the cloud mass 
of all disks is M c \ = 10 M Q . This is equivalent to a common 
column density of all clouds of 



_ 3 VLv tmh § 

cl ~4 ttGQ 



(23) 



A thick disk with -d /Q = 1/5 has clouds with column densities 
of N cl - 10 24 cm- 2 . 

5. Torus evolution 

The collisional disks described above can be identified with the 
observed 30 pc-scale gas concentrations observed by Davies et 
al. (2007) and Hicks et al. (2008) which arguably correspond 
to the large scale structure of AGN tori. 

It is now investigated how a collisional torus can evolve 
from a massive to a less massive state. To do so we assume that 
the gas mass of the torus is proportional to its mass accretion 
rate 



M gas - 2-3G-3tf3Q-3/-iM^ rotJ R = jM* 



(24) 



In this way we subsume all possible time dependencies of var- 
ious parameters and their correlations in the time evolution of 
M. Together with Eq.[l3]this yields 



ctR 2 







1 



64 jM c iG 2 M x 



(25) 



The expressions for the Toomre parameter, the turbulent veloc- 
ity, and the star formation rate are then 



i 



-f 2 G 2 M 2 fi M 2 



Wturb 



and 



M x ~ 



64 M ci j 2 G 2 n 
~d[R 2 



M 



2x 



(26) 



(27) 



(28) 



In this way the behavior of the transparency, Toomre parameter, 
thickness of the disk via H = Uturb / and star formation rate 
can be identified with the change of the external mass accretion 
rate. For our stationary equilibrium disks to be applicable we 
must require the changes of the external mass accretion to be 
slow, so that the whole disk can adjust to the changing external 
conditions. This time for adjustment t cq is approximately the 
ratio between the torus gas mass and mass accretion rate. For 



typical values of M gas ~ 10 7 M Q (Tab.|2|andM ~ lM Q yr _1 
(Tab.|4]i this leads to t cq ~ 10 Myr. This is smaller than the ob- 
served starburst ages (Tab.|2> for all AGNs except NGC 1097 
where the time of adjustment and the starburst age are compa- 
rable. 

5. 1. Evolution at constant thickness 

Whenever the disk thickness stays the same during its evolu- 
tion, t>turb =const. in time, and subsequently x — \ (Eq. |27l i, 
this implies that the external mass accretion rate stays at a high 
level. The gas mass, star formation rate, <&, and the Toomre Q 
parameter depend on the mass accretion rate M in the follow- 
ing way: 

M gas oc , M* oc M , i9 oc M"5 , Q oc . (29) 
The dependence of these parameters on the gas mass is then: 



gas 



(30) 



and the relation for the cloud mass and the volume filling factor 
are: 



M c i = const. , 4>v M„ 



(31) 



It is remarkable that the star formation rate is linearly coupled 
to the mass accretion rate in this scenario. 



5.2. Evolution at constant mass 

Alternatively, at later stages when the mass accretion rate 
is low and changes little, the torus mass may stay constant, 
Af gas =const, and subsequently x = (Eq. |27| >. The turbu- 
lent velocity dispersion, star formation rate, and the Toomre 
depend on the mass accretion rate M now in the following way: 



Uturb M 2 , M* = const. . fl = const. 



oc M' . (32) 



The relation for the cloud mass and the volume filling factor 
are: 



M ci = const. , 3>v cx M 2 



(33) 



At constant disk mass and decreasing external mass supply 
the disk will become geometrically thin without changing the 
transparency. 

An interesting result is that during both types of torus 
evolution — at constant turbulent velocity and at constant gas 
mass — the cloud mass does not change. 

6. Torus evolution scenarios 

In the picture of quasi-stationary equilibrium disks driven by 
cloud collisions the evolution will be determined by the ex- 
ternal mass accretion rate, i.e. the mass inflow from distances 
> 100 pc. We divide the torus evolution into three phases: 

- Phase I: Initial massive infall and formation of a turbulent, 
massive gas disk: 

An initial rapid infall of a large amount of gas , M gas ~ 
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10 7 M , within a short time (At < 1 Myr; M > 
10 M Q yr _1 ) leads to the formation of a massive (Q ~ 1), 
moderately thin (wturb /v TO t < 5) gas disk in which star for- 
mation proceeds. This phase will be recognized as a star- 
burst. The disk becomes turbulent and the turbulence is 
maintained by the energy input from stellar winds. After 
~ 10 Myr the first SN explode and will rapidly remove 
the initial intercloud medium from the disk. Only the most 
massive and densest clouds which are not Jeans unstable 
will survive. This leads in the following to a collisional 
torus. 

- Phase II: Torus evolution at constant turbulent velocity: 
During the first phase of its evolution the massive colli- 
sional torus stays thick. This implies that the mass accretion 
rate within the torus, and thus also the external mass accre- 
tion rate M, do not decrease significantly during this phase. 
The gas mass and the cloud collision rate (i co ii oc i? -1 ) de- 
crease with decreasing M, whereas Q increases with the 
square root of the external mass accretion rate. The star for- 
mation rate within the torus decreases with decreasing M 
(Eq.El. 

- Phase III: Torus evolution at constant gas mass: 

Once the external mass accretion rate has significantly de- 
creased, the torus evolves at constant gas mass. The veloc- 
ity dispersion and Q decrease with the square root of the 
external mass accretion rate, whereas the cloud collision 
rate (i co ii oc #) and the star formation rate stays constant 
(Eq. [32l i. The Circumnuclear Disk (CND) in the Galactic 
Center (Giisten et al. 1987) represents such a late stage 
(Q = 190, = 15; Vollmer et al. 2004) of the torus to 
disk evolution. 

Depending on the time evolution of the external mass accre- 
tion rate (from distances larger than 100 pc) we suggest three 
possible evolution scenarios (Fig.Q]i: 

- Scenario I: 

The torus never reaches Phase I, because (i) it is already 
clumpy from the very beginning, (ii) there is already star 
formation occurring at scales of ~ 10 pc from the galactic 
center; the associated SN explosions and/or winds remove 
the existing intercloud medium and/or inhibit the forma- 
tion of an intercloud medium. The two possibilities imply 
that the initial mass accretion rate is not very high. A third 
possibility (iii) is that the disk's velocity dispersion is pro- 
hibitively large to allow aQ = l disk. In this case star for- 
mation proceeds via cloud-cloud collisions. The star forma- 
tion rate is lower than that of a massive turbulent disk. The 
rate of momentum injection due to SNe and stellar winds is 

F = p*AA = 5 10 33 L. M * r ) dyne , (34) 

where p* is the pressure due to SNe and stellar winds ex- 
erted on a surface AA (Veilleux et al. 2005). Assuming the 
same initial turbulent velocity for the clouds and the in- 
tercloud medium, the intercloud medium can be removed 

from the disk if > PiM^v^turb' wnere ftM an d are 
the density and the volume filling factor of the intercloud 



medium. Assuming M» = 0.02 M Q yr _1 , B = 10 pc 
(Sect. [73]>, and v tulb = 50 kms -1 (Tab. |2) leads to 
niM&v 1 < 10 3 cm~ 3 . Since the disk volume averaged 
density of the clouds is n c /<I>v = 10 4 cm -3 (Tab. 0J, the 
intercloud space can only be cleared by SN explosions and 
stellar winds if the intercloud medium contains less than 
~ 10 % of the total disk mass. 

- Scenario II: 

Due to an initial, massive infall a massive (Q ~ 1) tur- 
bulent star-forming disk is formed (Phase I). The turbu- 
lence in this disk is maintained through the energy sup- 
ply by feedback from rapid star formation. The subsequent 
SN explosions destroy the disk structure after 10 Myr, i.e. 
the intercloud medium is removed leaving only the densest, 
most massive clouds which remain Jeans-stable. The disk 
becomes collisional and stays geometrically thick (Phase 
II). After ~ 100 Myr the mass accretion rate decreases and 
the disk becomes thin (Phase III) and ultimately transparent 
($ > 5). The time at which the torus changes from Phase II 
into Phase III depends on the time evolution of the external 
mass accretion rate in this scenario. 

- Scenario III: 

Due to an initial massive infall event, a massive (Q ~ 
1) turbulent star-forming disk appears (Phase I). As in 
Scenario II, SN explosions destroy the disk structure after 
~ 10 Myr, the intercloud medium is removed and only the 
densest, most massive clouds are left over, which are Jeans- 
stable. The disk becomes collisional and will stay thick as 
long as the external mass accretion rate is sufficiently high 
(Phase II). Due to a secondary massive and rapid gas infall 
a second massive (Q ~ 1) turbulent star-forming disk can 
form (Phase I), which evolves again into a collisional disk 
after ~ 10 Myr (Phase II). After ~ 100 Myr the mass ac- 
cretion rate has sufficiently decreased and the disk becomes 
thin (Phase III) and ultimately transparent ($ > 5). The 
time at which the torus changes from Phase II into Phase 
III depends again on the time evolution of the external mass 
accretion rate. 



7. Applying the observations 

The scenarios of Sect. [6] derived from our analytical models 
can now be compared with the VLT SINFONI observations of 
Davies et al. (2007) and Hicks et al. (2008). As has been shown 
in Sect. 14. II the present state of the disk can only be described 
consistently by a collisional disk model. The comparison with 
observations will allow us to derive the parameters of (i) the 
present torus, (ii) the initial torus immediately after Phase I, 
(iii) the massive (Q ~ 1) turbulent gas disk that gave rise to the 
initial starburst (scenario II and III), and (iii) the star formation 
efficiency of the collisional phase. 

The observables, i.e. the input parameters for our analytical 
model, are the radius from the galactic center R, the rotation 
velocity v ro t, the turbulent velocity dispersion v tul -h, the gas 
mass Mgas, the peak star formation rate during the initial star- 
burst Af* oak , and the age of the initial starburst tsB- We assume 
for all AGN, except NGC 1097, a starburst duration of 10 Myr, 
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Fig. 1. Schematic of torus evolution scenarios: Torus mass ac- 
cretion rate (M yr _1 ) is plotted as a function of time (Myr). 
Phase I: A massive (Q ~ 1) turbulent disk is associated with 
a starburst. Phase II: A collisional torus evolves at a constant 
turbulent velocity dispersion, i.e. thick torus evolution. Phase 
EH: The torus becomes thin and evolves at an approximately 
constant gas mass. 

i.e. the the starburst continues until the first SN explode. We 
only apply our model to a subsample of 6 nearby AGNs from 
Davies et al. (2007) for which these input parameters are suf- 
ficiently well known. The parameters for these objects can be 
found in Table |2] For each object, the table gives its Seyfert 
type and the radius within which the subsequent parameters 
apply. The dispersion Vt UI b is the value measured from the data 
(after accounting for instrumental broadening), while the rota- 
tion velocity v Iot is a Keplerian equivalent value. This means 
that it represents the rotation velocity that would be needed if 
ordered circular motions in a single plane supported the entire 
dynamical mass. It is therefore significantly greater than the 
measured rotation speed. The gas mass M gas is difficult to de- 
rive. Hicks et al. (2008) estimated it from a combination of di- 
agnostics, including the 1.3 mm CO luminosity, the 2.12 /1111H2 
1-0 S(l) luminosity, and a comparison to the gas-to-dynamical 
mass ratios in other spiral and starburst galaxies. The peak star 
formation rates M* e are derived from the starburst models 
used in Davies et al. (2007). They are simply the star formation 
rates required to form the young stars in a timescale of 10 Myr 
(in contrast to the time-averaged rates given in that paper). The 



Table 3. Equations to derive the disk/torus properties from ob- 
servations. Velocities are in units of kms -1 , radii in pc, masses 
in M Q , and star formation rates in M Q yr _1 . 



present collisional torus 



Q PI = (CrbMot) (M dyn /M| a r s ) 

tf pr = 1.9 x 10- 9 Q pr i?/(^ ot < r rb M cl ) 
M pr = 2« r rb ) 3 /(GQ pr tf pr ) 



massive turbulent gas disk 



?dis k = g674 ^ u ^ 



Mg£ k = 3.97 x W 5 !^, 11 R 
<S = 18.2M* 



initial collisional torus 



M init = M(M^f/Mgl s ) 2 



only exception is NGC 3783, for which we have adopted a 2- 
starburst model with ages of 1 10 Myr and 30 Myr. Such models 
were not considered by Davies et al. (2007) because of the lim- 
ited number of diagnostics. In our scenario with stronger theo- 
retical constraints on the models we need these additional de- 
grees of freedom to reproduce the observations. The final col- 
umn gives the age ige of the most recent starburst. All equa- 
tions to derive the properties of the present collisional torus, the 
massive (Q ~ 1) turbulent gas disk, and the initial collisional 
torus are given in Table [3] 

7. 1. The present collisional torus 

We recall the assumption that all gas clouds have a constant 
mass of M c \ = 10 M . Moreover the sound speed within the 
clouds is set to c s = 1.5 kmr 1 . The Toomre Q parameter 
of the present disks is directly calculated from Eq. [T] We then 
use the following expressions for d and the mass accretion rate 
M which follow from the expression derived for the turbulent 
velocity dispersion in Paper I w turb = {±G$Qf- 1 M)3 and 
Eq.Ql] 



= 7T7 c tG~ 1 Qtt- 1 M c 



(>4 



1 "turb 



and 



M = 2fvt ulh G- 1 Q- 1 ^' 1 



(35) 



(36) 
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name 


type 


R 


Viot 




M gas 










(PC) 


(kms" 1 ) 


(kms" 1 ) 


(10 6 M ) 


(Moyr- 1 ) 


(Myr) 


Circinus 


Sy2 


10 


93 


56 


2 


0.02 


80 


NGC 3783 


Syl 


30 


88 


35 


5 


0.3 


30 


NGC 3227 


Syl 


30 


173 


95 


20 


5.0 


40 


NGC 1068 


Sy2 


30 


184 


88 


20 


14.0 


200 


NGC 1097 


Syl 


30 


122 


59 


10 


1.0 


8 


NGC 7469 


Syl 


30 


123 


63 


10 


3.0 


110 



Table 4. Parameters of the present collisional tori. 



name 


M 

(Moyr" 1 ) 


Q 


•d 


$v 
10" 3 


Pel 

10 7 (cm -3 ) 


r c \ 

(PC) 


Circinus 


0.75 


5.7 


1.9 


1.2 


3.9 


0.01 


NGC 3783 


0.07 


4.1 


6.8 


0.2 


2.7 


0.01 


NGC 3227 


4.27 


5.5 


1.7 


0.5 


3.9 


0.01 


NGC 1068 


3.45 


5.4 


1.7 


1.2 


1.8 


0.02 


NGC 1097 


0.59 


4.8 


3.4 


0.4 


2.6 


0.02 


NGC 7469 


0.66 


5.1 


3.4 


0.6 


1.8 


0.02 



The torus mass accretion rates (Table @]i include the geomet- 
ric / = 0.1 factor introduced in Sec. 12.21 The resulting rates 
are still higher than the black hole mass accretion rates and 
the wind mass loss rates discussed in Sec. |SJ. The parameters 
derived in this way for our sample of 6 AGNs are shown in 
Table E] 

All tori show Toomre Q parameters between 4 and 6, 
i.e. they are massive thick tori. Moreover, 3 AGNs (Circinus, 
NGC 3227, and NGC 1068) have opaque tori (<Q < 2), 2 AGNs 
(NGC 1097 and NGC 7469) show moderately transparent tori 
(<# ~ 3 - 4), and 1 AGN (NGC 3783) has a transparent torus 
(i9 > 5). It is worth noting that two out of the 3 AGNs with 
opaque tori are classified as Sy2. We find the smallest mass ac- 
cretion rate for NGC 3783. The highest mass accretion rates 
(NGC 3227 and NGC 1068) are several ten times higher than 
that of NGC 3783. 

7.2. The massive (Q ~ 1) turbulent gas disk 

We now derive the parameters of the initial massive turbulent 
gas disk which gave rise to the initial starburst. For this we ap- 
ply the SN model (see Sect. |2|i where the energy source for 
maintaining turbulence in the disk are stellar winds. Their en- 
ergy input is comparable to that of SN explosions (MacLow & 
Klessen 2004). Therefore we do not need to change the formal- 
ism of the SN model. We further assume that the mass accre- 
tion rate equals the peak star formation rate given by Davies 
et al. (2007) (Table [2]) and <5 = 5, which is the Galactic value 
(Vollmer & Beckert 2003). 

We assume that for all AGNs, except for Circinus and 
NGC 3783, scenario II is valid. For Circinus we argue below 
that scenario I applies, because its peak star formation rate is 



a factor of more than 10 lower when compared to the other 
starbursts. This does not imply that star formation is not oc- 
curring, just that there was no initial massive accretion event. 
For NGC 3783 scenario III is more applicable, because a single 
starburst leads to an enormously high initial mass accretion rate 
compared to the present value, which we think is implausible. 
The double starburst we have adopted for NGC 3783, which 
is consistent with the observations of Davies et al. (2007), re- 
quires a first intense burst with M* =4.7 M Q yr~ 1 to have 
occurred 1 10 Myr ago, followed by a second distinct burst with 
M* = 0.3 MQyr" 1 only 30 Myr ago. 

The Toomre Q parameter, the total gas mass, and the tur- 
bulent velocity dispersion can then be calculated using the fol- 
lowing expressions: 

Q = 0.81 G-A M«r*a-*C*«rot , (37) 
M gas = 8.5 10- 5 G-^M *"<T"f , (38) 
and 

v tmh = 0.82GAm#<Ttt^ , (39) 

where £ = 4.6 10~ 8 (pc/yr) 2 (Vollmer & Beckert 2003). 

Using a peak star formation rate of M* = 0.02 M Q yr~ 1 
for Circinus leads to a gas mass of the massive turbulent disk 
which is smaller than that of the present collisional torus. We 
therefore conclude that scenario II (Fig. [TJ does not apply to 
this AGN. Instead, scenario I yields more appropriate results. 
The nuclear disk in Circinus did not experience a turbulent, 
supernovae and stellar wind driven Q ~ 1 disk. 

The gas masses of the turbulent starburst disks is between 
10% and 100% higher than the gas mass of the correspond- 
ing present collisional torus, except for Circinus. In the course 
of torus evolution (Fig. [TJ the loss of gas mass in the disk is 
moderate (up to a factor of 2; Tab. [2] and |5]l . The star forming 
disks are moderately thin (v ro t /^turb ~ R/H ~ 5-7) and their 
Toomre Q parameter is close to unity. 

7.3. The initial collisional torus 

Since all observed tori are thick (Table |2), they are all in Phase 
II of their evolution (see Sect. |6j, i.e. they evolve at constant 
thickness or velocity dispersion. This implies the following re- 
lation between the gas mass and the mass accretion rate: 

M oc Af 2 . (40) 
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Table 5. Parameters of the initial massive (Q ~ 1) turbulent 
gas disk that gave rise to the initial starburst. 



name 


M gas 

(10 6 M©) 


M 

(Moyr- 1 ) 


Utm-b 

(kms" 1 ) 


Q 


Circinus 


0.7 


0.02 


6 


1.9 


NGC 3783 


6.9 


0.3 


13 


1.1 


NGC 3227 


24.8 


5.0 


28 


1.3 


NGC 1068 


40.0 


14.0 


37 


1.2 


NGC 1097 


11.9 


1.0 


18 


1.2 


NGC 7469 


19.6 


3.0 


25 


1.0 



Table 6. Parameters of the initial collisional tori. 



1.0 












0.8 




01097 

\ s^GC3227 






- 


0.6 




>|«5C3783 








0.4 






^[rcinus 


5G7469 


^GCU 


0.2 
0.0 













50 100 150 200 

starburst age (Myr) 



name 


M 

(Moyr" 1 ) 


Q 





Circinus 


2.1 


3.4 


1.1 


NGC 3783 


0.1 


3.0 


5.0 


NGC 3227 


6.5 


4.4 


1.4 


NGC 1068 


13.5 


2.7 


0.9 


NGC 1097 


0.8 


4.0 


2.9 


NGC 7469 


2.6 


2.6 


1.7 



We can safely assume that the gas mass of the initial collisional 
torus at the end of Phase I is close to the gas mass of the massive 
(Q ~ 1) turbulent disk (Tabled. For Circinus we estimate the 
gas mass of the initial collisional torus by postulating that the 
disk was initially opaque, i.e. i? ~ 1. This leads to an initial gas 
mass of M gas = 3.4 10 6 M Q at the beginning of Phase II. 

For all other AGN a continuous transition from phase I to II 
provides the mass accretion rate of the initial collisional torus. 
The Q and $ parameters can then be calculated using the fol- 
lowing expressions (see Sect. |5j: 
-l 



Q oc M 



-l 

ras J 



■d oc M~ 



(41) 



Four of the initial collisional disks at the beginning of 
Phase II (Circinus, NGC 3227, NGC 1068, and NGC 7469) 
had Toomre Q parameters around Q — 3 and $ close to unity, 
i.e. they were massive and opaque. Two initial collisional disks 
(NGC 3783 and NGC 1097) had been massive and transparent. 



7.4. The evolution of the torus mass accretion rate 

In this section we investigate the evolution of the torus mass 
accretion rate with time in Phase II (see Sect. |6j. For this we 
plot in Fig. |2] the fraction between the mass accretion rate of 
the present collisional torus and that of the initial collisional 
torus (end of Phase I and beginning of Phase II) as a function 
of time for our sample. To determine the time that a given torus 
passed in Phase II, we place ourselves in scenario II (Fig.[TJi and 
adopt the starburst ages of Davies et al. (2007) for all AGNs 
except NGC 3783. We assume scenario III for this galaxy, i.e. 
the occurrence of two distinct starbursts. This leads to an es- 
timated age of the most recent starburst of 30 Myr. The solid 



Fig. 2. Ratio between the mass accretion rate of the present col- 
lisional torus and that of the initial collisional torus (end of 
Phase I and beginning of Phase II) as a function of time. The 
zero point in time is the end of Phase I ($ — 1 for Circinus). 

line is meant to guide the eye. Based on this plot we conclude 
that the mass accretion rate shows a slow monotonic decrease 
with time. It decreases to one fourth of its initial value in about 
150 Myr. 

7.5. The star formation efficiency of collisional tori 

Star formation in collisional disks is expected (see Sec. [3]) to be 
proportional to the cloud collision rate t~\, = Q/d: 



(42) 



where r\ is the star formation efficiency which we would like 
to determine. This is only possible for Circinus, because only 
in this case the peak star formation rate derived from obser- 
vations reflects the initial star formation rate of the collisional 
torus (scenario II). For all other AGNs the peak star formation 
rate is related to a massive (Q ~ 1) turbulent gas disk. During 
the torus evolution at constant thickness (Phase II) the star for- 
mation rate is proportional to the square of the gas mass 



M* oc Mz 



(43) 



The initial and present gas mass of the torus in Circinus are 
M™£ = 3.4 10 6 M Q and M gas = 2 10 6 M Q ; the peak star for- 
mation rate is M" ut — 0.02 MQyr -1 . This leads to an estimate 
of the present star formation rate of M* = 7 10~ 3 M Q yr _1 . 
The resulting star formation efficiency is 



V 



RM. 



init 



10" 



For a turbulent galactic disk the star formation law is 

A/* = ?7galMgas^ . 



(44) 



(45) 



Thus, one has to compare -q/d ~ 5 10~ 4 for a collisional disk 
with ?7 ga i = 0.017 derived for galactic SN driven turbulent gas 
disk. We conclude that the star formation in a collisional torus 
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is about 10 times less efficient than that of a turbulent gas disk 
of the same mass. In our sample of 6 AGN the star formation 
rate of the present collisional torus is small compared to the 
mass accretion rate. The average fraction between the mass ac- 
cretion rate and the star formation rate in these tori is 25 ± 7. 

7.6. Obscuring the starburst 

The stellar luminosities calculated by Davies et al. (2007; their 
Fig. 8) indicate that on scales of several tens of parsecs, the 
starbursts can only be weakly obscured. The reason for this 
is explained by Davies et al. (2006), and summarised here. 
If one assumes no obscuration, the starbust already comprises 
typically a few percent of the galaxy's bolometric luminosity. 
Even a moderate optical depth will lead to nearly all the opti- 
cal light being absorbed and re-radiated in the far-infrared; and 
in addition the scaling one derives for the starburst would in- 
crease. The limiting constraint is that the starburst luminosity 
cannot exceed the galaxy's bolometric luminosity. We estimate 
the resulting extinction (for screen models) to be in the range 
Ay — 4 to 8. This is consistent with a clumpy torus with area 
filling factors below unity (transparency z9 > 1) for radii larger 
than about 1 pc. The optical extinction is then mainly due to 
dust located in the foreground of the star forming region. In 
Sy2 galaxies, however, not only the central AGN is obscured 
by the torus, but also the inner part of the torus and embedded 
star forming regions are self-obscured by the extended torus. To 
verify this idea, we use the model of Beckert & Duschl (2004) 
and Honig et al. (2006) which is based on the formalism devel- 
oped by Vollmer et al. (2004) to estimate the extinction using 
(i) a screen model where the starburst occurred in the torus mid- 
plane behind half the torus, and (ii) a mixed model where the 
forming stars and the clouds share the same spatial distribution 
within the torus. As an example we use parameters most appro- 
priate for NGC 3227 or NGC 1068: v t urb/v Iot = 0.55, d = 2 
with an outer radius of 80 pc. Fig.[3]shows the mean number of 
clouds A along the line of sight through the torus for two torus 
inclinations in the case of a screen model. 



A 



Cud* 



(46) 



where Z co n = ■& H is the mean free path of the clouds and H the 
disk height. Obscuration with A > 3 leading to an likelihood 
of non-obscuration of less than e~ 3 occurs at a radii R < 20 pc 
for i = 60° and not at all for i — 40°. For lower inclinations 
the obscuration by clouds is ineffective. For i > 60° the ob- 
scuration pattern does not change much for geometrically thick 
tori. 

Fig. 2] shows the ratio between the extinction-free starburst 
emission and the emission with obscuration by intervening 
clouds in the mixed case 



A 



1 



(47) 



Complete absorption Iq / J ex t S> 1 is only reached in the inner- 
most part of the AGN torus (R <C 1 pc). 

The screen and the mixed models can be regarded as two 
extreme cases and the reality maybe somewhere in between 




X[pc] 




X[pc] 



Fig. 3. Projected mean number of clouds A along the line of 
sight through the torus for the screen model. Contour levels 
are (0.5,0.75,1,1.5,2,2.5,3,..). The thick contour is for A = 1 
corresponding to an non obscuration probability of e _1 . Upper 
panel: torus inclination i = 45°. Lower panel: torus inclination 
i = 60° (i measured from the torus axis). The size of the box 
is 100 pc. 



these two models. Most probably the starburst will be signifi- 
cantly obscured at galactocentric radii smaller than 1 pc. The 
extinction of star formation at the 10 pc scale by torus clouds 
is so low that additional extinction by extended dust lanes at 
larger radii (> 100 pc) in the galaxies is possible. We can there- 
fore conclude that our torus model is consistent with the small 
observed extinction of the central starburst. 
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Fig. 4. Ratio between the extinction-free starburst emission to 
the emission with cloud absorption Iq /Text for an torus incli- 
nation of i = 50°. Contour levels are (1.5,2,2.5). The thick 
contour is for Iq /I ex t = 2. It is the seconds outermost contour. 

8. Fueling the central engine 

The inner edge of the torus is thought to be set by the dust 
sublimation radius which is located at about 0.1-0.5 pc from 
the central black hole. In Sect. [7] we have derived the mass ac- 
cretion rate of the tori, i.e. the mass transport rate arriving at 
these inner edges. It is not clear what happens between the in- 
ner edge of the torus and the thin accretion disk around the 
black hole and what the relation is to the maser emission re- 
gions. These maser disks have sizes of about 0.1 pc. Between 
0.1 pc and 0.5 pc from the central black hole an X-ray heated 
wind is most likely formed (Krolik & Kriss 1995, 2001). Less 
dense and sheared clouds lose their dust by evaporation at the 
sublimation temperature are consequently ionized by the AGN 
X-ray emission, heated, and blown away in a line-driven wind. 
What remains of the dust or is newly formed in the wind to 
obscure the AGN at even higher inclination angles is undeter- 
mined in this scenario. The mass accretion rate onto the central 
black hole A/bh is thus the difference between the mass accre- 
tion within the torus M and the mass loss due to the AGN wind 

-A/wind^ 

Mbh =M- Mwi„d • (48) 

Typical wind mass loss rates are M W i„d = 0.03-0.3 M yr~ 1 
(Blustin et al. 2005, 2007). 

To investigate the relation between the mass accretion rate 
onto the central black hole and the present torus mass accretion 
rate, we plot their ratio as a function of the area filling factor 
$A = 4/3 The mass accretion rate onto the central black 
hole Mbh is derived from the AGN luminosity 

L = V Mbhc 2 , (49) 
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Fig. 5. Ratio between the mass accretion rate onto the central 
black hole and the present torus mass accretion rate as a func- 
tion the area filling factor, i.e. <1>a = 4/3/$. 

where the efficiency r) = 0.1 and c is the speed of light. The 
area filling factor is inversely proportional to the transparency 
of the torus (Fig. |5). Small area filling factors correspond to 
transparent tori. 

We observe a tentative trend in the sense that in opaque tori 
only 0.1-10 % of the torus mass accretion rate feeds the central 
black hole. On the other hand, if the torus is transparent, 10- 
100 % of the torus mass accretion rate is used to feed the central 
engine. 

The main differences between opaque (i? ~ 1) and trans- 
parent tori (d > 3) are the cloud collision rate and the cloud 
mean free path. Since there is no correlation with the collision 
rate t co \\ = $/f2, we suspect the mean free path to be respon- 
sible for the accretion efficiency of the tori. The different mean 
free paths between the clouds means that in transparent tori 
the AGN emission can reach all clouds of the torus, whereas 
in the opaque case it only reaches clouds at small galactocen- 
tric radii, i.e. close to the central source. Due to this radiation 
only the densest clouds can survive, the less denser clouds be- 
ing evaporated by a photodissiciation or X-ray dissociation re- 
gion (PDR/XDR). Vollmer & Duschl (2001) calculated the lo- 
cation of the ionization fronts in the cloud of the circumnuclear 
disk in the Galactic Center. These clouds have about the same 
mass as the AGN torus clouds (M c \ ~ 10 M Q ) and sizes of 
~ 0.1 pc. Whereas these clouds are close to the shear limit, 
the AGN torus clouds are ~ 100 times denser. Note however, 
that in the Galactic Center the source of ionization is a central 
cluster of about 40 O stars which is much weaker than an AGN. 

Vollmer & Duschl (2001) showed that the cloud radius due 
to the ionization front is given by 

r cl = 3.645 10 15 J^cf^d , (50) 

where Juv is the number of UV photons per cm 2 and sec, c\ 
is the sound speed in the ionized gas, and c s is the sound speed 

in the cloud. If one assumes that the cloud temperature is de- 

i i 

termined by the external radiation field c s ex Ta oc Jyy tne 



Vollmer, Beckert, & Davies: Starbursts and torus evolution in AGN 



13 



cloud radius does not change with the cloud's distance from the 
galactic center. 

We scale up the Galactic Center by a factor of 10 4 and 
obtain c\ = 7.6 kms -1 , c s — 3 kms _1 (R/(l pc)) _ 3, and 
J uv = 1.8 10 16 cm _2 s _1 (R/(l pc))~ 2 for a transparent torus. 
This leads to a cloud radius of r c \ = 0.03 pc, which is very 
close to the cloud radius obtained for our sample of AGN tori 
(Table |4]i. Thus, only dense, selfgravitating clouds can survive 
in an illuminated environment, i.e. a transparent torus. On the 
other hand, if the torus is opaque, less dense and larger clouds 
can survive. Since we derive only the properties of the dens- 
est and most massive clouds in our model which determine the 
physics of the outer torus, less dense and less massive clouds 
are still consistent with the model as long as they do not domi- 
nate the disk mass. As a consequence, the transparency values 
derived for the individual AGNs and summarized in Tableware 
upper limits. We therefore suggest that towards the inner edge 
of opaque tori clouds become less dense with densities close to 
the shear limit. These clouds, however, once they arrive at the 
inner edge of the torus, are (i) destroyed easily by the influence 
of shear and possible winds and are (ii) much easier ionized and 
evaporated by the AGN emission. In addition, for AGN close to 
the Eddington limit Honig & Beckert (2007) showed that dusty 
clouds experience a strong radiation pressure. For transparent 
tori only dense and compact clouds smaller than the shear limit 
will be found, while for opaque tori shadowing by clouds at 
the inner edge allows larger clouds to survive at intermediate 
distances. 

This tentative picture leads to two predictions for AGN tori 
that might be verified in the future: 

- the inner parts of opaque tori should harbor larger clouds 
and should therefore have larger volume filling factors, 

- the ratio of the torus mass accretion rate to the mass re- 
moved in winds generated at the inner boundary of tori 
should be lower for opaque tori. 



9. Uncertainties 

Up to this point we ignored any possible observational errors. 
In this section the influence of the quantities derived from ob- 
servations on the correlations shown in Fig. [2] and Fig. [5]is in- 
vestigated. We assume the following uncertainties: 

- For the gas mass (M gas ) - about a factor 2, because of the 
uncertainty in mass distribution and hence in interpreting 
the kinematics; 

- peak star formation rate (M* e ) - at least a factor 2-3 since 
the extinction is not well constrained; 

- starburst age (isB) - up to a factor 2, assuming that all diag- 
nostics are explained by a single burst, which we reject for 
NGC3783; 

- rotation velocity (v ro t) - ~ 40 %, since it is the equivalent 
Keplerian rotation velocity, and depends on both the mea- 
sured dispersion and inclination-corrected velocity. 

- turbulent velocity dispersion (wturb) - ~ 20 % assuming 
that the velocity dispersion is not strongly affected by a 
wind. 



All uncertainties of our derived quantities are thus dominated 
by the observational uncertainties on M gas and M*. For the 
torus mass accretion rate we have 



128 
M = —rfv. 



J "turb t 'rot^gas-^dyn C s 



(51) 



During the torus evolution the mass accretion rate is propor- 
tional to the square of the gas mass M oc M 2 as . The gas mass 
of the massive (Q ~ l) turbulent disk only depends on the peak 
star formation rate 



M, 



gas 



OC 



(MP eak ) 1 



(52) 



We took this gas mass as the initial gas mass of the collisional 
torus at the beginning of Phase II. Thus the ratio between the 
mass accretion rate of the present torus M and that of the ini- 
tial collisional torus M \& (Fig. O depends on the square of 
the present gas mass M gas and on the inverse of the peak star 
formation rate M| 



peak. 



M M, 
oc 



M 



old 



M, 



peak 



(53) 



On the other hand, the ratio between the mass accretion rate 
onto the black hole Mbh and the torus mass accretion rate M 
depend on the AGN luminosity L and on the inverse of the 
square of the gas mass M gas : 



Mi 



BH 



M 



3C 



M 2 

gas 



(54) 



An underestimation of the gas mass, which we assumed to be 
10 % of the dynamical mass, thus leads to a strong increase of 
the two ratios, while an underestimation of the peak star for- 
mation rate and the AGN luminosity leads to a decrease of the 
two ratios. In extreme cases both ratios can have errors up to a 
factor of 10. 

The starburst ages have uncertainties of a factor of 2. The 
area filling factor is proportional to the inverse of the gas mass 
3>v oc M^g (Eq.[35l). Thus, the associated uncertainty is also a 
factor of 2. 

Despite these relatively large uncertainties we believe that 
the correlations shown in Fig. [2] and [5] are real. All system- 
atic errors, like an overestimate of the gas mass, would alter 
but not destroy the correlations, as long as the errors for dif- 
ferent targets are not random. However, to make our findings 
more robust more spectroscopic VLT SINFONI observations 
with high spatial resolution of AGNs are needed. Most impor- 
tantly, future ALMA high resolution CO line observations are 
necessary to determine the total gas mass of the tori with an 
uncertainty of ~ 10 — 30 %. 

10. A holistic view of torus evolution in AGN 

VLT SINFONI observations of the close environments (~ 
30 pc) of a sample of nearby AGNs by Davies et al. (2007) 
showed that thick gas tori and recent central starbursts with 
ages smaller than 100 Myr are ubiquitous. We compare differ- 
ent clumpy accretion disk models to these observations: 
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- fully gravitational turbulent disks where the turbulence is 
maintained by the energy input from the gravitational po- 
tential via mass accretion, 

- supernova and stellar wind driven turbulent disks (SN 
model) where the turbulence is maintained by stellar winds 
and supernova explosions, 

- collisional disks where the orbital motion is randomized by 
partially elastic collisions which also allow mass transport 
to the center and angular momentum redistribution. 

Whereas the measured rotation velocity, turbulent velocity dis- 
persion, and gas mass tell us about the current state of the gas 
torus, the peak star formation rate and the age of the starburst 
provide information on the past appearance of the torus. We 
assume that the physical properties of the torus are mainly de- 
termined by external mass accretion from scales of ~ 100 pc. 

The result of this work is a time sequence for the torus evo- 
lution. Present tori appear to be collisional and geometrically 
thick whereas the tori giving rise to the starburst in the past 
are of turbulent nature and relatively thin. The torus evolution 
can be divided into 3 phases depending on the external mass 
accretion rate: 

- Phase I: initial massive infall: formation of a massive tur- 
bulent stellar wind-driven gas disk with Q <~ 1, 

- Phase II: decreasing, but still high mass accretion rate: col- 
lisional thick torus, 

- Phase III: decreasing, now low mass accretion rate: colli- 
sional thin torus. 

Phase I is short (~ 10 Myr). Once the SN explode, they remove 
the intercloud medium and clear the disk leaving behind only 
the densest clouds. The result is a collisional torus. All tori dis- 
cussed in this paper are interpreted to be in phase II. Therefore 
this phase can last for more than 100 Myr. The transparency 
of these tori depend on their Toomre parameter, the turbulent 
velocity, and the internal pressure of the clouds (temperature 
and magnetic fields). If there is a second short massive infall, 
the torus can again switch into phase I. We suggest that this 
has happened in the case of NGC 3783. On the other hand, 
an initial massive gas infall is not mandatory as the case of 
Circinus illustrates. Once the external mass accretion rate has 
significantly decreased, the torus becomes thin and transparent. 
The Circumnuclear Disk in the Galactic Center represents this 
phase. 

In addition we conclude that 

1 . The massive turbulent stellar wind-driven gas disk (phase 
I) gives rise to a starburst with a star formation rate that 
equals the disk mass accretion rate. 

2. The collisional torus also forms stars, but with an efficiency 
which is about 10 % that of the turbulent disk. 

3. During the evolution of the thick collisional torus (phase II) 
the decrease of the mass accretion rate is slow (a factor of 
4 within 150 Myr). 

4. The present collisional tori (phase II) do not significantly 
obscure the stellar populations born during the recent star- 
burst. This is the case even for Sy2 nuclei. 

5. We find a tentative correlation between the area filling fac- 
tor of the clouds and the ratio between the mass accretion 



rate onto the central black hole and the torus mass accre- 
tion rate. AGNs with a high area filling factor lose more 
than 90 % of the mass transported through the torus in a 
thermal torus wind at the inner edge. We suggest that this is 
due to shadowing of the gas clouds from the central AGN 
throughout most the extended torus which allows the exis- 
tence of low density clouds whose densities are close to the 
shear limit. 

Future VLT SINFONI and ALMA observations will be 
necessary to confirm our proposed torus evolution sequence. 
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